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Abstract. A technique to simulate the spex-mixer/mill system as a macroscopic ensemble rather than a pure dynamical sys- 
tem is proposed. The treatment is suitable especially for comminution processes generating the nanomaterial up to nanometers 
scale where the geometrical displacements are unobservable. It is argued that the method is simple and experimentally verifi- 
able since relevant physical observables can be extracted using partition function without solving the equation of motions. 
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INTRODUCTION 

The comminution processes in recent years attract the attention among scientists and engineers due to the increasing 
demand of ultrafine powders for nano technology applications in many areas. The demand then requires the improve- 
ment of comminution equipments like ball mills, roller mills and so on. Unfortunately, the development of such 
comminution equipments always contains a lot of uncertainties due to a wide range of unknown parameters. These, in 
fact, lead to significant statistical errors. In order to overcome such problems, several models have been developed to 
quantitatively describe comminution process in various types of mills [1, 2, 3, 4, 5]. 

On the other hand, mathematical modeling and simulation may provide prior information and constraint to the 
unknown parameter ranges which should be useful to develop more optimized experimental strategy in comminution 
processes. However, in most cases of mathematical models, the physical observables like grain-size etc are extracted 
from a set of equation of motions (EOM). Such EOM's are considered to govern as complete as possible the 
dynamics of the system, from the mechanical motions to the evolution of grain-size distribution. This approach is 
obviously suffered from the nonlinearities of the equations under consideration, and then the requirement of high 
computational power to solve them numerically. This fact often discourages a quantitative and deterministic approach 
for the simulation of such system. These nonlinear effects like chaotic behavior of the sphere motions within the mill 
encourages some works modeling the system using semi-empirical approaches [6]. However most of semi-empirical 
models require a large number of experimental data based on prior observations [7], or measured variables obtained 
from simulation results by other authors [8, 9]. 

More empirical approach is based on the physically realistic modelization of the ball mill system [10, 11]. This 
approach in general deals with three aspects : (1) evaluation of milling bodies dynamics and energetic inputs transferred 
to powders; (2) description of the effects of such inputs on powders breakage; (3) description of powders evolution in 
terms of particle size distribution [12]. In a recent work [12], a comprehensive study on this line for the case of spex 
mixer / mill was performed by deploying the 3D simulation for milling bodies motion, and the population balance 
method to describe the granulometric evolution. Then, both are related through the energetic inputs in the population 
balance formalism which is obtained from the simulation of milling bodies motion. 

In this paper we propose a novel model and approach combining the deterministic approach for milling bodies 
motion, and the statistical approach to relate them with considerable macroscopic physical observables. This work is 
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devoted to overcome the following problems in conventional approach : 



• Experimentally it is almost impossible to trace the geometrical displacements of all matters in a vial with proper 
time resolution to verify the models which are based on the classical EOM. This problem is getting worse as one 
simulates a system with matters at few nanometers scale with comparable size of time-space resolution. 

• Solving a set of EOMs numerically, and then performing a simulation with high accuracy (enough time resolution) 
require huge efforts on both computing capacity and running time. 

• Taking into account the external circumstances around the vial like electromagnetic field and so forth. This might 
be interesting when one considers a comminution process which can reach the level of few nanometers. 

Therefore this work is intended to provide a tool for a nanometer system in a vial by developing direct relations 
between the vial internal dynamics with some external physical observables which should be easier to measure. We 
should remark here that the vial internal dynamics is yet described empirically using physical modelization approach. 

Further, rather solving the EOM's governing the whole dynamics, we use the hamiltonian approach to model all 
relevant interactions, and extract the physical observables through partition function by considering the system as a 
canonical ensemble in finite temperature. As a consequence, instead of observing the geometrical evolution of matters 
in term of time in a ball mill, we can investigate the particle number distribution in term of temperature. This introduces 
a novel method relating the internal dynamics with the macroscopic physical parameter like temperature, rather than 
time and geometrical displacements which are in most cases difficult to realize. An example of numerical simulation 
is given for the case of ball mill with a structure similar to the well-known spex mixer / mill. 



THE MODEL 

The whole system is modeled empirically using hamiltonian method. First we construct the total hamiltonian describ- 
ing the dynamics in the ball mills. It is further followed by formulating the partition function and extracting the relevant 
thermodynamics observables. 



The dynamics 

In our model, the dynamics of each 'matter' in the system, i.e. balls and powders inside the vial, is described by a 
hamiltonian H m (r,t). The index m denotes the powder (p) or ball (b) and r — (x,y,z). The hamiltonian contains some 
terms representing all relevant interactions working on the matters inside the system as follow, 

H m = Ho + Vm-m + Vm-v + Vm-m' + ^ext , ( 1 ) 

with v denotes the vial, while Hq is the free matter hamiltonian, that is the kinetic term, 

1 «m 

^^Eiouf, (2) 

where n m is the matter number, m m and p m are the matter mass and momentum. Throughout the paper we assume that 
the mass or size evolution of matters is uniform for the same matters. 

The matter self-interaction V m _ m , the matter-vial interaction V m _ v and the interaction between different matters 
may be induced by, for instance, Coulomb (y Coul ) and impact (y lm P) potentials, 

i=lj=l J >J 
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(4) 



with Q m is the matter charge, while m,m' : v,p,b and n is the unit normal vector. These potentials are considered 
describing the mechanical and static electrical properties of the matters. The first potentials should in fact represent 



the whole classical dynamics among the matters, i.e. the impact forces among balls and powders. This form will be 
clarified soon below. The Coulomb potential disappears if the interacting matters have neutral charges. Also it works 
only in a short range of distance. Therefore, it should be negligible for nanometers scale of physics as in our case. The 
impact forces between the vial surface and balls or powders are treated in the same way using Eq. (3) by considering 
that the surface is constructed from a set of fixed spheres [12]. This is inline with the simulation in the last section 
where we discretize the vial volume in small spheres with a comparable size as the desired powder size, i.e. few tens 
nanometers at the largest. 

On the other hand, instead of Eq. (3) we can consider a simpler 'effective potential' like the harmonic oscillator 
^m-m'^) = 2 ^mm' Ar 2 to represent the whole mechanical properties in terms of 'effective coupling' k mm i. In this case, 
if m m » m m i, the potential can be rewritten as V m _ m / = ^ m m 0) m 2 Ar 2 since CDm = \/ k mm i / m m . Actually this is the case 
of ball-powder interaction since m\, ^> ra p by the order of namely (9(10 2 ). Nevertheless, in contrast with its simplicity, 
y osc absorbs the time dependency and also interesting physical parameters characterizing the material properties like 
viscoelasticity, Young modulus etc. The time dependency is important to directly relate the system temperature with 
physical observables through finite temperature partition function as shown in the next subsection. Therefore, in this 
paper we take the impact potential to represent the mechanical properties in the system. 

The impact potential in Eq. (3) is induced by the impact force consists of normal and tangential components [12], 
Fj(r,() =r ] ,(f,()|F | J i] ,(r,(). The normal component may be written [13], 
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(5) 



Here the first term is the elastic part based on the Hertz contact law, and the second term is the dissipative one that 
takes into account material viscosity. Y mm / is the Young modulus and v mm / represents the Poisson ratio of the sphere 
material. The term , — (/? m /? m /)/(/? m + /? m /) represents the effective radius, while <i; mm / = R m +R m , — \r m — r m >\ is 
the displacement with ^ m is the radius of interacting matter. A is a dissipative parameter [13, 14, 15], 
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The viscous constants r] m and rj m i relate the dissipative stress tensor to the deformation tensor [13, 14]. 

There are several proposed formulations for the tangential component F^ m ,(r,t). However it always assumes that 

the material slide upon each other in the case where the Coulomb condition jj, F^ m , < F^ m , holds, otherwise 

some viscous frictions occur [16]. In particular it follows that F^ m ,(r,t) °< m mm'> wri ere the effective mass is m^ m , = 
m m m m i I (m m + m m i) [12]. Obviously, in our case with large mass discrepancy between powder and ball, the tangential 
impact force may be neglected for a good approximation. So, let us from now consider the normal component 
dominated impact force, that is F^,(f,t) ~ F^ m ,(r,t). This result simply yields the impact potential as written in 
Eq. (3) due to the Euler-Lagrange equation, 
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(7) 



since the dependency on matter velocity appears only in the tangential component ,(r,t ) [12]. 

Beside the interactions among the matters itself, it is also possible to take into account the external potentials 
working on the whole system. For instance, in the dynamics of ball, the gravitational potential, 



(8) 
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might be important, with G is the gravitational constant. On the other hand, this should be less important for the 
powder dynamics due to its tiny size. However, we may take into account the effect of external electromagnetic field 
surrounding the system to the charged matter. The potential is induced by the Lorentz force, F™ = Q m (E + v m x B), 
which leads to, 
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and satisfies Eq. (7). <j> and A are the electromagnetic scalar and vector potentials related to the electric and magnetic 
fields by E = —V<p — dA/dt and B = V x A. The inclusion of electromagnetic potential shifts the kinetic term in Eq. 
(2) as follow, 
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From now, let us focus only on the dynamics of powder which is our main interest in the sense of comminution 
process. From Eqs. (1), (2), (3), (4) and (9), the total hamiltonian for the powder in our model is, 
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for 6p 7^ 0- The last two potentials represent the total impact potential among powders; powders and vial; powders and 
balls respectively. Obviously we do not need to take into account the ball self-interaction V™^ nor ball-vial interaction 
y™P -p n j s j s ac tually the advantage of using hamiltonian method. 



Physical observables 

As mentioned briefly in introduction, the greatest advantage of deploying the hamiltonian method is one can extract 
some physical observables without solving the EOM's governing the system. Instead, in a canocical ensemble we 
consider the partition function of the model, 



Yldpidr, exp - dtH n 



(12) 



where j3 = l/(ksT) with kg and T are the Boltzman constant and absolute temperature. Having partition function at 
hand, we can obtain some thermodynamics quantities in the system through relations, 
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P m = — 



(14) 



for free energy and, 

dF m _ _Fm 
dV ~ V ' 
for pressure in a vial with volume V. 

In order to see the contributions from the interactions, it is more convenient to consider the normalized partition 
function, 
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and further, 
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(16) 



Performing the integral over time (t ), we immediately obtain a temperature dependent partition function, while the 
integrals over /?, are the decoupled gaussian integral which can be easily calculated. In the case of Eq. (1 1) it gives, 



p/2 r «1 



(17) 



and, 



where the interaction hamiltonian is, 
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Obviously, only the scalar potential of external electromagnetic field contributes to the total energy of system under 
consideration. In other words, we can conclude here that in our model the magnetic field B does not influence the ball 
mill system, but the electric field E does. 

Moreover, we can perform the integration over time (f ) and <^, 7 to obtain further, 
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From this result, the thermodynamics observables are clearly not affected with the dissipative term, i.e. the second 
term in Eq. (5). After performing same integration we obtain, 
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(21) 



respectively with, 



J (=1 L 13 "P,yj=l i-V PP V 1 -" ! Pm:b,vi=lj=l 1 "V 

Eq. (2 1 ) provides a general behavior for temperature-dependent pressure in the model, while the geometrical structure 
and motion of vial is absorbed in the function & . From Eq. (21) clearly the physically meaningful regions are for 
< T < {2mp%kB)~ l and T > T^. The later is equivalent to the condition, 



& <k B Tt tl ln(2m p nk B T (tl ) , 



(23) 



and T t h is always greater than (2m p Kkg) 1 . 



SUMMARY 

The simulation can in principle be performed to deal with any types of ball mills. This can be accomplished by 
replacing the coordinate system of ball mill under consideration represented by & in Eq. (22) to the appropriate one 
which represents its geometrical motion. 

We have proposed and discussed a novel model and approach for comminution processes of nanomaterial using ball 
mill equipments. The study is focused on investigating the hamiltonian for a ball mill system, and relating it to relevant 
physical observables through partition function. 

Detail simulation and its numerical results are still under progress and will be reported elsewhere. 
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